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1.0 


Introduce ion 


Satellite infrared (IR) radiometric measurements pro- 
vide a wealth of data on atmospheric parameters like temperature 
and humidity profiles. A major problem in this technique is the 
presence of clouds in the field of view of the radiometer. SMITH 
et al (1970) have used a statistical method for detecting cloud 
contamination. The same technique with some modifications is being 
used as part of an automated set of procedures to produce global 
maps of sea surface temperature on an operational basis by NOAA. 
When both IR and visible pictues of the same scene are available, 
the visible data can be used (SHENK and S0L0M0NS0N, 1972) to 
pick out cloudy pixels. A possible problem in this method is the 
mis-registration of the IR and visible images. Figure 1 shows 
how such a mis-registration can show that a visible pixel is clear 
while the corresponding IR pixel is cloudy. 

In the present study the effect of mis-registration on 
cloud brightness threshold is investigated by simulating radio- 
metric data as observed from a spin stablized synchronous 
satellite such as the SMS. Clouds are introduced randomly and 
a bi-directional reflectance model is used to create radiance 
data from clouds and ocean. As part of this study a theoretical 
and an empirical reflectance model are compared. 

2 . 0 Brightness Normalization 

2 . 1 Introductio n 

Brightness normalization is the technique of adjustment 
of visible radiance data measured under different illumination 
and viewing conditions (Figure 2) to a standard configuration of 
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Registration offset between visible (continuous 
lines) and IR (dashed lines) . Irregular contour 
represents cloud. 



- 2 - 


ORthtkal PAGE IS 
* OUALITY 





« » 




Figure 2 Geometry of Reflection 
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solar illumination anti sensor viewing. This is necessary in 
applications where intercomparison of cloud albedo has to be 
made. The simplest normalization procedure would be to assume 
that clouds are perfectly diffuse reflectors which obey Lambert's 
law. However, studies of the angular distribution of reflected 
solar radiation (RUFF et al, 1970, BRENNAN and BANDEEN, 1970, 

RASCHKE et al, 1973) have shown that the reflectance of clouds 
depends on the relative geometry of both the sun and the satellite 
with respect to the cloud. In such a case a bidirectional re- 
flectance p has to be defined (SIKULA and VONOERHAAR) as follows: 

Reflected radiation from 
the direction 0, 4/ , 

pU, 0, <P) - (1) 

Incident radiation from 
the direction £. 

where £, 0. <l> are defined in Figure 1. It can be shown that the 
directional reflectance r (£) is related to p by the relation 

r(C) - itp (0 (2) 

At present two anisotropic models, one empirical 
(SIKULA) and VONDI RHAAR, 1972) and the other theoretical (MOSHER, 
197“) are available for normalizing cloud brightness. 

2. 2 E mpirical Model 

SIKULA and VONDERHAAR (1972) have synthesized an empirical 
bidirectional model by using observations of cloud reflectances by 
BRENNAN and BANDEEN (1970) as well as available summaries presented 
by RASCHKE et al (1971). Figure 3 shows a vertical cross-section 
of the model. 
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The empirical model has the following defects: 

a. It has a restricted data base In the sense that re- 
duction of variance in reflectance will require a 
much larger number of observations than have been 
actually used. 

b. For purposes of cloud height estimation the optical 
thickness will not be available as a parameter of 
the model. 

It also has the following advantages: 

a. It is simple to use. 

b. It can be updated continually, as and when new ob- 
servations become available. 

c. Its empirical nature frees it from errors due to 
assumptions inherent in the theoretical model. In 
other words, it is 'true to life'. 


2 . 3 The oretical Mode l 

The light reflected from a cloud depends on the following 
variables (MOSHER, 1974) : 



a. Shape and size distribute i n of particles. 

b. Number density of scattering particles. 

c. Cloud thickness 

d. Angular configuration of sun, cloud, and satellite 

e. Shape of the cloud 
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MOSHER (1974) has developed a model of cloud reflectances 
by theoretically calculating the multiple scattering of light In 
a cloud for a range of optical thicknesses of the cloud. The 
multiple scattering program is th* one developed by HANSEN (1971). 
MOSHER (1974) has computed the intensity of reflected light from 
the cloud for a range of elevation angles of the sun and the 
satellite as well as the relative azimuth angle between them. The 
model is thus four dimensional. Three angles and the optical 
thickness of the cloud are the parameters of the model. Given 
the satellite-cloud-s»in geometry and the measured radiance, it 
will then be possible to pick out the optical thickness from the 
model . 

The theoretical model has the following defects: 

a. The theoretical model neglects both Rayleigh scat- 
tering and molecular absorption. At small optical 
thicknesses, the former may have to be taken into 
account. There is a possibility (CURRAN, 1975) that 
absorption due to ozone may be large enough to be 
taken into account. 

b. The multiple scattering computation using the 
doubling method cannot take the finite horizontal 
dimension of the cloud into account. There is evi- 
dence from Monte Carlo calculations (MCKEE and COX, 

1974, DAVIES, 1975, MCKEE and COX, 1975) that there 
is a considerable leak of radiation through the 
sides of the cloud. The theoretical model corrects 
for this effect by an arbitrary adjustment of a 


scaling factor used to convert input digital counts 
into radiance data. The factor was adjusted to give 
'reasonable results' when the emissivlty of cirrus 
clouds was conputcd. 

c. Finally, the method of computing percentage cloud 

cover is open to question. However, the vexing pro- 
blem of partial cloud cover in a resolution element 
i 8 one which is common to all satellite radiometric 
data. 

The model has the following advantages: 

a. It is also simple to use. 

b. It yields optical thickness as an extra parameter 
which makes it more suitable for cloud height esti- 
mation. 


2. A Conclusions 

The theoretical approach becomes 'empirical' when the 
calibration factor is adjusted. The validity of such an adjustment 
has to be tested by calibrating the visible data independently. 

A good check on the MOSHER model would be e comparison with results 
obtained from Monte Carlo calculation ror the same angular con- 
figuration for finite clouds. The calibration factor for very 
thick convective clouds of finite horizontal dimensions can then 
be computed from the Monte Carlo results. 

When these checks show the theoretical model to be a 
valid one, it will be obviously superior to the empirical model. 


♦ 
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3.0 Simulation of Visible Radiances as Observed from a 

Synchronous Satellite 

Simulation of radiance data Involves the following 
procedures : 

1. The starting and ending numbers of lines and pixels 
of the grid under consideration are assumed. The 
dimensions of the grid were taken to be 64 by 64. 

2. In the case of a spin stabulized synchronous satellite 
like the SMS the mapping raster is formed by a combina- 
tion of satellite spin motion and a step action of the 
scanner optics. Assuming the orbit and attitude para- 
meters of satellite, to be the nominal values of SMS-l, 
the rotation matrix for conversion from the satellite 
coordinates (Figure 4) to earth-fixed coordinates 
(Figure 5) is computed. 

3. For each pixel the pointing vector from the spacecraft 
to the pixel is computed in satellite coordinates, and 
then transformed to earth-fixed coordinates through the 
rotation matrix. The piercepoint of this vector with 
the earth ellipsoid is then computed. The latitude and 
longitude of each pixel are then easily obtained. 

4. Cloud cover was simulated by assuming that one quarter 
of the 64 x 64 pixels had clouds with randomly varying 
areas in them. Uniformly distributed random numbers 
varying in magnitude from 0 to 1 were generated for 
this purpose. 

5. The bidirectional reflectance pattern of SIKULA and 
VONDERHAAR (1972) was used to compute the albedo A 
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observed by the satellite where A is given by 
A - ^ P -— + °‘ 7, ffy( L " C ^ COS? 


The 

numbers .03 and .7 refer to the 

zenith albedo of 

the 

ocean and 

clouds respectively. 

P and Q are defined 

by 





P - fo 

<C> 



r 

o 

TFff) 

(4) 

and 

Q “ Ic 

(O 

(5) 


r 

c 




Where r is the isotropic reflectance and the subscripts 
C and 0 refer to the cloud and ocean respectively. P is 
the anisotropic reflectance and XX and YY are the aniso- 
tropic factors defined by 


XX - r c (t> 

(6) 

*o c (? ,0 .CO 


YY = r o 

(7) 

"P 0 (CO ,|H 



P, Q, XX and YY have been tabluated by SIKULA and VON 
DER HAAR (1972). 

The digital count corresponding tc the albedo sensed 
by the radiometer was taken as D where 

D = 62.8 Y~A (8) 

which is the nominal calibration function for the VISSR 
instrument (BRISTOR, 1975). Observational noise was added to 
the digital counts by generating Gaussian random numbers wirh 
a standard deviation of 1 count. A sample of the simulated 
digital counts over a clear ocean area can be seen in 
Figure 6. 
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FIGURE 6 

DIGITAL COUNTS OF VISIBLE RADIANCE 
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4.0 Determination of Visible Brightness Threshold for Clouds 

o 

The digital counts of albedo over the 64 x 64 grid are 
normalized by dividing by cosf, , under the assumption of isotropic 
reflectance. This is because the anisotropic reflectance for 
ocean is different from that for clouds. Unless the percentage 
cloudiness of a pixel is known, it is not possible to use the 
anisotropic ref lectance model for normalization. A histogram of 
the normalized reflectance was then created. Figure 7 shows a 
computer plot of the histogram. The horizontal scale shows both 
the sequential interval class and the normalized reflectance 
corresponding to the interval. Interval 1 contains values less 
than 8 and interval 20 all valves exceeding 20. The occurrence 
frequency in each interval is given in the top line. From Figure 8 
it may be seen that the histogram is asymmetric. The peak in the 
histogram corresponds to the predominant clear ocean area. Points 
to the left of the peak are due to the assumed observational noise. 
SHENK and SAL0M0NS0N (1972) determined the width of the despersion 
Ar due to noise on ocean reflectances from the histogram. The 
brightness threshold was then set as r^FAr where r^ is the modal 
value of reflectance. 
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5 . 0 Error duo to Isotropic Normalization 

As pointed out in Section A.O, it is not possible to use 
anisotropic reflectance models for normalizing the brightness of 
partially cloudy pixels. It is therefore necessary to use the 
isotropic reflectance assumption in normalizing the brightness of 
visible pixels. In order to study the possible error due to this 
assumption, a visible image matrix of 64 x 64 pixels was generated 
with the empirical anisotropic reflectance model as described in 
Section 3.0. The visible radiance counts were normalized with 
the assumption of isotropic reflectance. The histogram of 
normalized reflectances is shown in Figure 7. The visible thres- 
hold value (Section 4.0) can be seen to be 11. This threshold 
value was used to identify 'clear' pixels. These 'clear' pixels 
were then tested against the original simulated matrix to determine 
how many of them were really clear. The number of cloudy pixels 
erroneously classified as clear and the total cloudy area wrongly 
classified were thus identified. This procedure was repeated for 
successively larger thresholds. Table I shows the results corres- 
ponding to a grid simulated almost vertically below the satellite 
when the sun’s zenith angle was about 20°. It can be seen from 
Table I that corresponding to a threshold of 11 only 16 pixels out 
of a total of 4096 pixels were erroneously classified as cloudy. 

The total area of cloudiness in these 16 pixels was about a seventh 
of a pixel. These 16 pixels that were classified as clear contained 
clouds that were very small in extent and as such were comparable 
to observational noise. 

Obviously, increasing the visible threshold will increase 
the number of cloudy pixels wrongly classified as clear. However, 
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for a random cloud model the visible threshold can be increased 


beyond the r + Ar value (Section 4.0) without incurring serious 
error on a statistical basis. From Table I it can be seen that 
increasing the threshold from 11 to 18 will increase the total 
cloudy area wrongly classified as clear to 4 pixels out of a total 
of 4096 pixels. 

TABLE I 


C - 19.6 Q - 58.4 f- 60 . 6 


Visible 

Threshold 

(Normalized 

Radiance 

Counts) 

Number of 
Cloudy Pixels 
Erroneously 
Classified 
as Clear 

Total Area 

Classified Erroneously 
as Clear 

(Units of Pixels) 

11 

16 

.14 

12 

22 

.29 

13 

35 

.72 

14 

45 

1.19 

15 

53 

1.57 

16 

61 

2.02 

17 

73 

2.90 

18 

86 

3.96 

19 

94 

4.72 

20 

107 

6.11 


However, if there was a very thin layer of cloud which 
increased the brightness by say 7 counts, the error in the surface 
temperature estimation due to the relaxed threshold will be considerable. 

To investigate the effect of the geometric configuration 
of satellite and sun with respect to the area under consideration, 
a 64 x 64 matrix was generated nearer the limb, were the values of 
£, 0 and <j / differed considerably from the previous case. The 
histogram of normalized reflectances for this case is shown in 
Figure 8. The visible brightness threshold in this case can be 
seen to be 18. The number of cloudy pixels erroneously classified 
as clear corresponding to various thresholds is shown in Table II. 
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In this case also the error due to relaxation of the threshold can 
be seen to be small. Although the Isotropic normalization wLll 
lead to different values of the threshold for different geometric 
configurations of sun and satellite, the use of the brightness 
threshold value does not introduce significant error. This is 
because the brightness threshold selection technique is a statistical 
one and error due to isotropic normalization acts like a constant 
bias over small areas. 


TABI.E II 


C - 58.9 

0 - 48.9 

V - 8.3 

Visible 

Number of 

Total Area 

Threshold 

W xels 

Erroneous ly 

(Normalized 

Erroneously 

Classified 

Radiance 

Classified 

as Clear 

Counts) 

as Clear 

(Units of Pixels) 

18 

8 

.06 

19 

8 

.06 

20 

12 

.12 

21 

12 

.12 

22 

14 

.18 

23 

21 

.31 

24 

27 

.55 

25 

35 

.89 

26 

39 

1.10 

27 

47 

1.41 


6 . 0 Effects of Spatial Misregistrat i on Between Visible and 

IR Data . 

In order to use the brightness threshold of clouds 
obtained from the visible data to identify clear column IR pixels, 
it is necessary for the IR and visible pixels matrices to be in 
exact spatial alignment. If the IR and visible images are mis- 
aligned then the identification of clear column IR pixels from the 
visible data loads to error. 

The effect of misregistration on the cloud filtering 
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process was tested by creating a visible radiance matrix 
dimensions 64 x 64 pixels with a constant offset along each line 
from the IR matrix. In other words, the simulated misregistration 
between the IR and visible matrices was onlv along the lines and 
none across them. Such a misregistration has been observed the 
SMS VISSR data. 

For a given offset the number of cloudy pixels wrongly 
classified as clear was determined for different visible thresholds. 
Table III gives the results for an offset of one pixel along each 
line. It may be seen from Table III that there is a sharp increase 
in the area wrongly classified as clear. The error is clearly 
unacceptable. To illustrate this effect further, the offset 
between the IR and visible matrices was increased to 4 pixels along 
the line only. The errors resulting from this can be seen in 
Table IV. It is obvious from Table IV that the increase in 
registration offset increases the error enormously. 



TABLE III 


Visible 


Number of 


Total Area 


Threshold 
(Norma lized 
Radiance 


Clouds 
Erroneous ly 
Classif ied 


Erroneous ly 
Classified as 
Clear (In 


Count) as Clear Units of Pixels) 


11 

72 

41.2 

12 

78 

44.3 

13 

90 

49.2 

14 

99 

54.5 

15 

105 

58.2 

16 

113 

61.8 

17 

124 

67.0 

18 

136 

73.1 


It has to be remembered that the above results refer to 
a cloud model that is random. In real data cloud strata are clearly 
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Identifiable and only scattered clouds within the field of view 
cannot be visually identified. A random cloud model can represent 
one realization of such scattered clouds within the field of view. 
Hie results shown in Table III and IV are valid for such cases. 


TABLE IV 



Visible 

Threshold 

(Normalized 

Radiance 


Number 
of Clouds 
Erroneously 
Classified as 


Total Area 
Erroneous ly 
Classified as 
Clear 


Count) Clear 


11 

239 

125.5 

12 

245 

129.2 

13 

256 

133.7 

14 

264 

140. 1 

13 

269 

143.6 

16 

277 

147.0 

17 

287 

151.1 

18 

295 

154 . 3 


The conclusion to be drawn from these results is that 
the IR and visible data have to be spatially registered to an 
accuracy of at least one pixel for the purpose of cloud filtering 
of IR data by using a visible brightness threshold. 
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